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For arbitrary networks of random masses connected by random springs, we define a general 
strong disorder real-space renormalization (RG) approach that generalizes the procedures introduced 
previously by Hastings [Phys. Rev. Lett. 90, 148702 (2003)] and by Amir, Oreg and fmry [Phys. 
Rev. Lett. 105, 070601 (2010)] respectively. The principle is to eliminate iteratively the elementary 
oscillating mode of highest frequency associated with either a mass or a spring constant. To explain 
the accuracy of the strong disorder RG rules, we compare with the Aoki RG rules that are exact at 
fixed frequency. 

I. INTRODUCTION 

In the field of disordered systems, the problem of random masses connected by random springs is very old since it 
has been introduced by Dyson [l| even before the classical Anderson localization paper Q. After studies concerning 
the one-dimensional case (see the review Q and references therein), an analysis of disordered elastic media via a 
non-linear sigma-model [|| has predicted results similar to the scaling theory of Anderson localization Q : all finite- 
frequency phonons are localized in dimension d < 2, whereas there exists a finite critical frequency transition between 
delocalized and localized modes in dimension d > 2. More recent discussions on the similarities and differences with 
Anderson localization of electrons can be found in 041] and references therein. 

Besides the case of regular lattices, one may also consider the case of more general networks with various physical 
motivations. Two types of strong disorder real-space renormalization (RG) procedures have been previously intro- 
duced: Hastings Q has proposed to eliminate iteratively masses, whereas Amir, Oreg and Imry (Io| have proposed 
to eliminate iteratively springs. In these two studies, the results have been found to be extremely accurate at low 
energies. In the present paper, we unify these two particular procedures into a single framework, where both masses 
or springs can be iteratively eliminated in a consistent way. In addition, we explain the accuracy of strong disorder 
RG procedures at low frequency in these models, by a comparison with the Aoki RG rules that are exact at fixed 
fre quen cy (this type of renormalization has been first introduced by Aoki in the context of Anderson localization in 

[Dint). 

The paper is organized as follows. Section |TT] presents the models and notations. In section IIII1 we describe the 
strong disorder RG rules consisting in the iterative elimination of the site or the link associated with the highest 
frequency, and explain the relations with the previous approaches 0, [Io| . In section IIV1 we compare with the Aoki 
RG rules that are exact at fixed frequency. Our conclusions are summarized in section IVl 

II. DISORDERED NETWORK OF MASSES CONNECTED BY SPRINGS 

We consider an arbitrary network where sites are indexed by i. With each site is associated a random mass m^, 
and with each link < i.j > is associated a random spring constant Ki j — Kj t i. 

The scalar phonon model is defined by the following harmonic Hamiltonian for the scalar displacements Ui(t) 

ff = Ey«?+E %"< ",r (l) 

i <i,j> 

The scalar assumption is very standard to simplify the analysis Q and means physically that vibrations along different 
directions are decoupled. Equivalently, the model can be defined by the equations of motion 



rriiUi 



= ^ Kjj (uj - Uj) (2) 



As stressed in 0, [T3] , this scalar phonon formulation can be used to study various other relevant physical models : 
(i) the case m, = 1 with random Kjj can be used to study the Laplacian on the network with disordered couplings 
between nodes [si, [13]. For instance in [lfjj|. the sites were drawn at random in a d-dimensional cube and the Kij was 
exponentially decaying in the Euclidean distance r^j. In [9[, the motivation was coming from the field of complex 
networks (see the recent reviews [HI])- 
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(ii) the case rrii = 1/rii where rij represents the number of nodes connected to i and Kij = 1 when are 
neighbors on the network (with Kij — otherwise) corresponds to the usual random walk on the network. Here the 
'disorder' comes only from the network heterogeneous structure. 

III. STRONG DISORDER RENORMALIZATION (R.G.) APPROACH 

Strong disorder renormalization (see fill ] for a review) is a very specific type of RG that has been first developed in 
the field of quantum spins : the RG rules of Ma and Dasgupta [16| have been put on a firm ground by D.S. Fisher who 
introduced the crucial idea of "infinite disorder" fixed point where the method becomes asymptotically exact, and 
who computed explicitly exact critical exponents and scaling functions for one-dimensional disordered quantum spin 
chains [13] • This method has thus generated a lot of activity for various disordered quantum models : see the review 
[lj| for works before 2004, as well as more recent developments concerning entanglemen t fl8,[l9ll . superfluid-insulator 
transition of disordered bosons [2(| , Bose- Einstein condensation [2l[ , dissipation effects 22l. 



23|, Hubbard model [24 



anyonic chains [25j], fractal lattices l26jl . studies on the link with extremal statistics [27|, and implementation of new 
very efficient numerical procedures |28| . Strong disorder renormalization has been also successfully applied to various 
classical disordered dynamical models, such as random walks in random media either in one dimension [2^, [3(j, on 
strips [31], in two dimensions [Hj], or in the presence of absorbers [33|, reaction-diffusion (34|, coarsening dynamics of 
classical s pin chains [35[ , trap models (36l , random vibrational networks [t| [l(| , contact processes [H, [3i| , zero range 
processes [401 ] and exclusion processes |41| , non-equilibrium dynamics of polymers or interfaces in random media [42j j , 
statistics of valleys in configuration space of disordered systems [43| , oscillator synchronisation [lij , and extreme value 
statistics of various stochastic processes |45( . 

In this section, we describe a strong disorder renormalization procedure for the network of masses connected by 
springs of Eq. [1] and [5] that generalizes the previous procedures proposed in d,[l(|. The principle of this RG procedure 
is to eliminate iteratively the highest frequency present in the system. The elementary oscillating modes are of two 
types as we now discuss. 

A. Oscillating modes associated with a single spring 

Let us first consider the problem of a single spring between two masses m\ and m 2 , when all other spring constants 
are supposed to be negligeable. The equations of motion 

mi?ii + Kifl (ui - u 2 ) =0 (3) 
m 2 it2 + Ki t 2 (u 2 ~~ Ui) = (4) 

can be analyzed by introducing the center-of-mass displacement uq and the relative displacement u re i 

m\U\ + m 2 U2 ,_ s 

u G = ■ (5) 

mi + m2 

Urel =U2—U\ (6) 

as well as their corresponding masses 

m G = mi + m 2 (7) 
TO1TO2 

mi + 7712 



m re i = ■ (8) 



The equations of motion of Eq. [4] are then decoupled 

m G ug =0 (9) 

m re l Urel + K 1,2 Urel = (10) 

One obtains the zero- mode corresponding to the global motion of the center of mass, and the oscillating model of the 
relative coordinate at the frequency Q re i given by 

n?, 2 = ^ = K1.2 (— + —) (ii) 



m re i \nii TO2 
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If this frequency is high, this mode will not be excited by the low-energy modes of the whole system : one may 
then eliminate the corresponding relative mode, i.e. we set 

Urel = U 2 — Ul = (12) 

The masses mi and are then merged into their center of mass G of mass tuq — ra\ + m 2 and of displacement 
uq = ui = U2- The equation of motion for the center of mass now takes into account the forces of other points 
i (1,2) that were connected either to (1) or to (2) or to both 

m G u G = ^ ( K i,i + K i?) ( u i ~ u g) (13) 

#(1,2) 

that had been neglected in Eq. [10] in the single spring analysis. The equation of motion for i ^ (1,2) becomes 

miili = ^ K i,n (u n ~ Ui) + (K it i + K it2 ) (ug ~ Ui) (14) 

71^(1,2) 

In summary, the elimination of a spring K\i associated with a high frequency fii^ (see Eq. II 1 p consists in the 
merging of the two sites 1 and 2 into a single site representing their center of mass G of mass tog = Tn% + mi- The 
approximation of Eq. [T^] yields that the springs constant connected to G are given by to the following renormalization 
rule 

K^ = (K i>1+ K it2 ) (15) 

whereas the masses m, are unchanged. 



B. Oscillating mode associated with a single mass 



The problem of a single mass mo in the environment of the other masses rrij that are supposed to be fixed ( 
Ui{t) = Ui independent of time) is of course very simple : from the equation of motion 



m u 



Uq 



one obtains that uo(t) will oscillate around the equilibrium value 

eq _ Sj K 0,i u i 



with the frequency S7q given by 



= — E Ko,i 



(16) 



(17) 



(18) 



If this frequency is high, this mode will not be excited in the low-energy modes of the whole system : one may then 
eliminate it. This amounts to say that «o(t) will adiabatically follow the slow motion of its neighbors 



uf m \t) 



Y,i K ,iUi(t) 

En K 0,n 



(19) 



The elimination of the mass mo leads to the following renormalized of the equation of motion for the other masses 



rriiUi 



£ 



K 



(uj - m) 



(20) 
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In summary, the elimination of a mass mo associated with a high frequency Hq (see Eq. [T8|) leads via the adiabatic 
approximation of Eq. [T5]to the following renormalizations for the spring constants between two neighbors of the 
mass mo 



whereas the masses m, are unchanged. 



C. Statement of the strong disorder renormalization procedure 



The above analysis suggest the following strong disorder renormalization procedure : 

(i) at a given RG step, there are a certain number of masses m^ and a certain number of spring constants Ki j = Kji 
(with i =/= j). 

With each spring Kij, one associates the frequency Qjj given by 



1 1 

1 

ra, rrij 

With each mass m^, one associates the frequency £li given by 



n 2 itj = K itj ( — • — ) (22) 



The renormalization scale is defined as the highest frequency f2 remaining in the system 

Q = max[{Qi,Q itj }} (24) 

among all frequencies associated with masses or spring constants. 

(ii) Decimation of the mode associated with the highest frequency f2 : 

(ii-a) If the highest frequency D, — Qi J0 is associated with the spring constant Ki j : the two masses mj and 
mj are replaced by their center of mass G(io,jo) of mass 

m G (i ,jo) = m *o + m j ( 25 ) 
and the spring constants linked to iq or to jo or to both are replaced by spring constants linked to their center of mass 

K 3,G(io,jo) = K 3,i + K J,h ( 26 ) 

(ii-b) If the highest frequency — Qi is associated with the mass «o : the mass rrii is eliminated, and the spring 
constants between two neighbors of (?o) are renormalized according to 

K K 

K i,j =K ^ + x^l? — ( 27 > 

(iii) Return to point (i) to update the frequencies associated with the surviving renormalized mass and to the 
surviving renormalized springs. 



D. Relation with the RG rules proposed by Amir, Oreg and Imry |l0j 

In [loj . Amir, Oreg and Imry have proposed to choose at each step the largest spring constant K% j and then to 
apply RG rules that coincide with Eqs [25] and [26j They insist that one should choose the largest spring constant 
independently of the masses, and not the highest corresponding frequency that we introduced in Eq. 1221 We do not 
agree on this point, since we believe that the frequency scale is the only well founded physical variable to define a 
consistent strong disorder procedure for phonons. We believe that for the specific application studied in [Toj . the 
distribution of spring constants is so broad whereas all masses are initially equal, that the two choices are probably 
equivalent. However for other applications, we believe that the frequency choice is the relevant one to compare the 
possible excitations of various sub-systems. At the beginning of the RG procedure, the considered sub-systems are 
made of a single mass or of a pair of masses connected by a spring, but as the RG proceeds, the considered sub-systems 
contain more and more initial masses. The frequency criterion allows to construct in a consistent way the appropriate 
renormalized modes that will respond to an exterior low-energy excitation. 
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E. Relation with the RG rules proposed by Hastings [9] 



In Q , Hastings first performs the similarity transformation [l[ 

Vi 



to obtain a symmetric operator from the equation of motion of Eq. [2J 



(28) 



with 



3 

Li , = giL. (30) 

The first RG procedure proposed in consists in choosing the highest Lij at each step : since Li : i = Qf of Eq. 
|2"31 this is equivalent to choose the highest frequency associated with masses only (and the frequencies fif^ of Eq. [22] 
associated with links are not considered for this choice). The RG rule corresponding to the elimination of Li Qt i reads 
Q (after correction of the typo concerning the sign of the second term in 



L^T = L U - L f Ll °' J (31) 



which is equivalent to the RG rule of Eq. [27J for the spring constants and to keeping the masses unchanged for the 
remaining sites. 

The second RG procedure proposed in Q still consists in choosing the highest L^j at each step, but to treat 
separately its neighbor j with the highest \Lij\ : the RG then consists in the diagonalization of this two-by-two 
submatrix and to project out the eigenvector of highest eigenvalue. In @, the advantages of this procedure with 
respect to the previous rule of Eq. [31] were to allow the increase of the connectivity and the variations of the masses. 
However these properties are also present in the simpler rules of Eqs 1251 and 1261 concerning the decimation of a link. 
We thus believe that the RG rules described in section IIII CI that consider both the decimations of sites and links 
is physically clearer and is able to renormalize at the same time spring constants, masses and connectivity of the 
network in a consistent way. 



F. Validity of the strong disorder renormalization procedure 

This renormalization procedure will be consistent if the renormalization scale f2 defined in point (i) decreases at each 
step, i.e. if the new generated frequencies are smaller than the decimated frequency f2. Moreover from what happens 
in other fields [15| . one expects that the procedure will become asymptotically exact at small f2, if the renormalized 
distribution of the frequencies become broader and broader as —¥ 0. For an arbitrary network with an arbitrary 
disorder of masses and spring constants, it is not easy to know a priori if the procedure will remain consistent, and if 
the procedure will become asymptotically exact at low frequency. However, this can be checked numerically for each 
case of interest. It turns out that in the various cases studied numerically in [t| [l^ the results have been found to 
be extremely accurate at low energies. In the following section, we explain that this accuracy could come from the 
coincidence with the zero-frequency limit of Aoki exact RG rules. 



IV. COMPARISON WITH AOKI EXACT RENORMALIZATION PROCEDURE AT FIXED 

FREQUENCY 

For Anderson localization models, Aoki pdl4l3| has introduced an exact real-space renormalization procedure at 
fixed energy which preserves the Green functions of the remaining sites. This procedure has been further studied for 
one-particle models in [46M48j | . It has been also extended in configuration space for two-particle models [49| and for 
manybody localization models [50j . This method can be extended to any eigenvalue equation, so it can be applied to 
the scalar phonon problem at fixed frequency as we now explain. 
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A. Equations of motion at fixed frequency 

Since the equations of motion of Eq. [2] are linear, one may also analyse the dynamics in terms of oscillating modes 
u u of a given fixed frequency ui : 

~ m l w 2 u w (z) = J~)Kj,j (u u (j) - & u (i)) ( 32 ) 

3 



B. Aoki renormalization procedure 
Eq. [32] for an arbitrary site i = io can be used to eliminate Uw(io) with 



Then the equations of the remaining sites i ^ io remain of the same form as Eq. 



with the renormalized frequency-dependent springs constants 



i \u) = Kijiu) + — — — 2 TT (.35) 



and the renormalized frequency-dependent masses 

rtii (u) = rrii{uj) + 5 p: (36) 

Let us stress the two essential differences with the strong disorder RG rules discussed in the previous section : 

(i) Aoki RG procedure is exact for any order in the choice of the decimated sites io, whereas in strong disorder RG 
procedure it is essential to choose at each step the mode corresponding to the highest elementary frequency. 

(ii) The price to pay is that Aoki RG concerns a fixed frequency Li, and that the renormalized parameters contain 
this frequency. So for each frequency, one should restart at the very beginning the RG procedure to compute the new 
renormalized couplings. On the contrary, the renormalized masses and springs of the strong disorder RG procedure 
do not contain the frequency. 

So these two RG procedures are completely different in spirit and in practice. But it turns out that they become 
similar if one consider the Aoki RG rules in the limit of zero-frequency. 



C. Limit of zero- frequency of Aoki RG rules 

In the limit of zero frequency u — > 0, the elimination of io corresponds to the RG rule (Eq. I35|) 

(0) = K id (0) + ^0(0)^(0) (37) 



which coincides with the strong disorder RG rule of Eq. [57] Note that the other RG rule of Eq. [M] concerning the 
masses disappear in the limit lu — > since all masses appear always together with ui 2 in the equation of motion (Eq. 
EH). ' 

This coincidence between the strong disorder RG rules and the zero-frequency Aoki exact RG rules may explain 
the numerical accuracy of the strong disorder RG rules at low frequency found in [i| IToj . Moreover, since Aoki 
procedure is exact for any order of the decimated sites, one expects that 'bad decimations' during the strong disorder 
RG ('bad decimations' meaning decimations when the eliminated frequency is not well separated from the others) 
will actually be corrected in the later stages of the renormalization when the neighbors will themselves be decimated. 
This phenomenon has already been found in other applications of strong disorder RG, in particular for quantum spin 
chains (see Appendix E of [17j |) and for the construction of valleys in configuration space for disordered systems (see 
section 3.6 of Si). 
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V. CONCLUSIONS AND PERSPECTIVES 



In summary, we have proposed a strong disorder renormalization procedure for the random phonon problem on 
arbitrary networks, that generalizes previous approaches where only site decimations @ or only link decimations [Toj 
had been considered. We believe that the iterative elimination of the mode of highest frequency associated with either 
a mass or a spring constant is the appropriate formulation to study cases where both processes can a priori occur 
during the RG process (in analogy with quantum spin chains in random transverse fields where one decimates either 
sites or bonds with an energy gap criterion (TR Il7|). To explain the accuracy of the RG rules found in the various 
examples studied numerically in 0, [Io| , we have compared with the Aoki RG rules that are exact at fixed frequency. 

More generally, strong disorder RG rules are naturally defined on arbitrary networks, since even if the starting 
point is a regular lattice of dimension d > 1, the RG rules changes the connectivity and thus destroys the initial 
regular structure (the one-dimensional lattice is usually the only case where the renormalized structure remains one- 
dimensional). So they can be applied to disordered systems defined on complex networks. However it turns out that 
complex networks constitute by themselves a disordered structure, so that even non-disordered models defined on 
non-regular networks can present Griffiths phases [5l[ or localization effects [52] . It would be thus interesting in the 
future to study whether strong disorder RG rules could be appropriate to study non-disordered systems on complex 
networks that are sufficiently inhomogeneous. It turns out that among the various RG procedures that have been 
recently introduced for complex networks |53r56[, a recent proposal [57| consists in a purely sequential algorithm 
where at each step, a node is selected at random to perform an elementary RG transformation : however since the 
nodes are not all equivalent, a natural question is whether the selection of a node at random is really the good 
choice, or whether one should choose the node or the link according to a maximal criterion, as in strong disorder RG 
procedures. 
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